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Abstract 

Up to now the raise and peel model was the single known example of a one- 
dimensional stochastic process where one can observe conformal invariance. The model 
has one-parameter. Depending on its value one has a gapped phase, a critical point 
where one has conformal invariance and a gapless phase with changing values of the 
dynamical critical exponent z. In this model, adsorption is local but desorption is 
not. The raise and strip model presented here in which desorption is also nonlocal, 
has the same phase diagram. The critical exponents are different as are some physical 
properties of the model. Our study suggest the possible existence of a whole class of 
stochastic models in which one can observe conformal invariance. 

1 Introduction 

There is a long list of papers on one- dimensional interface growth models (see [1] and [2] 
for reviews). In most of them the stochastic processes are local. Two typical gapless phases 
are encountered: the Edwards- Wilkinson phase [3] where the dynamical critical exponent 
z = 2 and the Kardar-Parisi-Zhang phase |1] where 2 = 3/2. If in a model one has an 
Edwards- Wilkinson phase by introducing an asymmetry one can get a Kardar-Parisi-Zhang 
phase. By introducing "friction" (some processes oppose growth) one gets a gapped phase 
(see the flip-flop model described below). The interplay of asymmetry and "friction" is nicely 
displayed in [5]. 

As far as we know the first example of a nonlocal growth model is the Derrida and 
Vannimenus' study of the interface in weakly disordered systems [6]. Much later another 
model with nonlocal rates appeared in the literature, the raise and peel model (RPM) [7]. 
This is a one-parameter (denoted by u) dependent model which also displays three situations, 
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analogously to local models. For small values of u one has a gapped phase, for n = 1 one has 
a gapless phase with z = 1 (not z = 2) and for n > 1 one has a gapless phase with varying 
values of z (this corresponds to the KPZ phase in local models). Dyck (special RSOS) 
paths describe the interface. Adsorption is local but desorption is nonlocal. The desorption 
processes look a bit artificial, they come from the algebraic background of the model at 
u = 1. What is special in this model is that for u = 1 one can do analytic calculations and 
show that the model is conformal invariant. Moreover, the stationary state of the model has 
fascination combinatorial properties [HlHlllO]. Recently, another model (the peak adjusted 
raise and peel model (PARPM)) was introduced [11]. The adsorption and desorption process 
are like in the RPM but the rates depend on the number of peaks in the Dyck paths. This 
makes the rates dependent on the size of the system. A new parameter p was introduced 
such that if p = 1, one recovers the RPM at m = 1. It was shown that conformal invariance is 
maintained in whole domain of p. In this paper we present a new model, the raise and strip 
model (RSM) which is again a one-parameter dependent model, with local adsorption and 
nonlocal desorption processes, the rules for the latter being much simpler. The configuration 
space is the same as in the RPM. Our aim was to see if the main ingredient to get conformal 
invariance is the existence of nonlocal processes of a special kind. We were aware that the 
price to pay is loosing integrability and that one had to use Monte Carlo simulations on large 
lattices to get results. We have indeed observed that the RSM has a phase diagram similar 
to the one seen in the RPM. At the conformal invariant point, the critical exponents are 
different. Since other models having the same structure as the RSM can easily be defined, it 
is plausible to assume that there is a whole class of models with conformal invariance which 
should be studied. 

The paper is organized as follows. In Section 2 we defined the observables for models 
defined on Dyck paths. These observables are used in the description of the properties of 
the models. 

The fiip-fiop, the raise and peel and raise and strip models are defined in Section 3. The 
flip-flop model is the local version of the RSM, it was studied in order to see the effect of 
nonlocality introduced in the RSM. 

The flip-flop model is presented in Section 4. One shows, using known results from 
combinatorics and Monte Carlo simulations that one has a critical point with a dynamic 
critical exponent z = 2 which separates two gapped phases. 

The raise and peel model is shortly reviewed in Section 5. 

The main results of our research are given in Section 6 in which we present not only 
the properties of the stationary states of the RSM but also the time dependent phenomena. 
An interesting new phenomenon occurs if the parameter u is larger the Uc- The system 
stays gapless with varying values of z (like in the RPM) but unlike the RPM where in the 
stationary state, the average height increases logarithmically with the size of the system, in 
the RSM, the average height profile is a triangle with a height of the order of the system 
size. 

For completeness, a variant of the flip-flop model in which the configuration space is 
changed, is presented in the Appendix. 

Our conclusions can be found in Section 7 
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Figure 1: An example of a Dyck path for L = 14. There are four contact points, three 
clusters, four peaks and three valleys. The substrate profile is shown in blue. 

2 Observables for models defined on Dyck paths 

We consider an open one-dimensional system with L + 1 sites {L even). A Dyck path is a 
special restricted solid-on-solid (RSOS) configuration defined as follows. We attach to each 
site i integer heights hi which obey RSOS rules: 

h,+i-hi = ±l, (z = 0,l,...,L-l), (2.1) 

with the constraints: 

ho = hL = 0, (2.2) 
/i, >0, (^ = 0,1,...,L). (2.3) 

There are 

Z{L) = L\/{L/2)\{L/2 + l)\ (2.4) 

configurations of this kind. 

A Dyck path can be seen as an interface separating a film of tilted tiles deposited on a 
substrate, from a rarefied gas of tiles (see Fig. [1]). The substrate corresponds to the special 
Dyck path defined as h2k = 0, /i2fc+i = 1 (/c = 0, 1, . . . , L/2 — 1). 

If hj = at the site j, one has a contact point (there are four contact points in Fig. [1]). 
Between two consecutive contact points one has a cluster (there are three clusters in Fig. [1]). 
The slope at the site i is s{i) = (/ij+i — /ii_i)/2. If Sj = {hi > hi-i) one has a peak. If 
Si = {hi < hi-i) one has a valley (there are four peaks and three valleys in Fig. [T]). 

It is useful to use the known mapping between an RSOS path and the configuration space 
of a one-dimensional hopping model with exclusion. To an upwards step hi — > one 
associates a particle, to an downwards step hi — hi^i < one associates a vacancy. The 
constraints f l2.2p -( !273l) are translated into two conditions on the particles- vacancies configu- 
rations. The first constraint (12. 2 p implies that the number of particles is equal to the number 
of vacancies. The second constraint (12. 3p is obviously nonlocal: the number of particles on 
the left side of any bond has to be larger or equal to the number of vacancies on the left 
side of the same bond. For example, the configuration XXOXOO is acceptable but not 
OOXOXX {X is a particle, O is a vacancy). The Dyck path shown in Fig. [T] corresponds 
to the configuration XXXOOOXOXXOXOO. 
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Figure 2: The flip-flop model. When a tile hits the interface, the following processes occur: 
The tiles c and d are reflected, tile h triggers a local desorption and the tile a which hits a 
valley, is adsorbed. 

The raise and fall stochastic models to be described below give the probabilities of the 
various Dyck paths and one is interested in average values of observables. 

One obvious observable is h{i, L, t) which is the average height at the site i for a system 
of size L and time t. The average density of contact points which is a function of the site i, 
the size L of the system and t, will be denoted by g{i, L,t). The average density of clusters 
equal to the average number of clusters divided by L, will be denoted by K{L, t). When one 
considers the stationary states of the models, the time dependence will be dropped in the 
notation. For example, g{i,L,t) will become g{i,L). It turns out that the average density 
of peaks and valleys T{L,t), equal to the average total number of peaks and valleys divided 
by L, plays an important role in studying the properties of the models. 

3 Raise and fall models 

We present three stochastic models defined in the configuration space of Dyck paths. In one 
of the models (the flip-flop model) the processes are local, in the other two (the raise and 
peel and raise and strip models), the adsorption processes are nonlocal. As we are going to 
see it is the nonlocality of the rates which is relevant in getting new physics. 

The models depend on one parameter u which is the ratio of adsorption and desorption 
rates. One uses sequential updating. At each time step, with a probability 1/(L — 1) a tile 
hits the Dyck path at a site i = 1, 2, . . . , L — 1. The effects of the hits are different in the 
three models. 

The flip-flop model (FFM) is an extension of the freely jointed chain model of a random 
coil polymer [T2]. In this model both the adsorption and desorption processes are local. The 
nonlocality comes only from the constraint (12. 3 p which defines the Dyck paths and, unlike 
the two other models, not from the rates. The model is defined by the following rules: 

If a tile hits a valley (sj = and hi < /ij_i), with a rate u it sticks to the site and the 
valley becomes a peak (see tile a in Fig. [2]). If a tile hits a peak (sj = and hi > hi^i > 0), 
with a rate equal to one, the tile at the peak gets desorbed and the peak becomes a valley 
(see tile b in Fig. [2]). If the tile hits a site i and Sj 7^ 0, the tile is reflected with no changes 
in the profile (see tiles c ans d in Fig. [2]). In this model the rates are local. In the raise 
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Figure 3: The raise and peel model. When tile a hits the interface it triggers a nonlocal 
desorption. The desorption can be local only if the tile hits the substrate (tile d). If the tile 
hits a peak, it is reflected (tile b). If a tile hits a valley, like in the flip-flop model, the tile is 
adsorbed (tile c). 

and strip model to be described below the desorption rates are nonlocal and it is interesting 
therefore to compare the models. 

The raise and peel model (RPM) was intensively studied [7]. We present it here in order 
to clarify the effects of nonlocality on the physics of the models. 

The main merit of the RPM is that for n = 1, it is integrable and conformal invariant. 
The essential features of the model are local adsorption and nonlocal desorption processes of 
a particular kind. The rules for the latter come from the algebraic structure (the Temperley- 
Lieb algebra) behind the model. These rules do not easily generalize for other configuration 
spaces where we expect to be able, based on algebraic considerations (using Hecke algebras), 
to define stochastic models which are also conformal invariant. An example is the case of 
restricted Motzkin paths [13]. The RPM is defined by the following rules: 

Depending on the slope Sj = (/li+i — /ij_i)/2 at the site i, the following processes can 
occur: 

1) Sj = and hi < (tile c in Fig. [3]). The tile hits a local minimum and with a rate u is 
adsorbed (/ij — > /ij + 2). 

2) Sj = and hi > hi^i (tile h in Fig [3]). The tile hits a peak and is reflected. 

3) Sj = 1 (tile a in fig. [3]). With a rate one the tile is reflected after triggering the desorption 
{hj — hj — 2) of a layer of 6 — 1 tiles from the segment {j = z + 1, . . . , -i + 6 — 1} where 
hj > hi = /ij+f,. 

4) Si = —1 (tile d in Fig. [3]). With a rate one, the tile is reflected after triggering the 
desorption [hj — )■ hj — 2) of a layer of 6 — 1 tiles belonging to the segment {j = i — b + 
1, . . . ,i — 1} where hj > hi = hi_fy. 

The raise and strip model (RSM) was conceived in order to keep the main features of the 
RPM (local adsorption and nonlocal desorption) using simpler rules and to see if one can 
recover conformal invariance. The price to pay is of course, lack of integrability and other 
magic properties of the RPM. Keeping in mind that in a Dyck path each valley is followed 
by a peak and that a peak is surrounded by two consecutive valleys, the RSM is defined by 
the following rules: 

If a tile hits a valley (sj = and hi < hi_i), with a rate u it sticks to the site and the 
valley becomes a peak (see tile a of Fig. Hj). If a tile hits a peak (sj = and hi > > 0), 
surrounded by two valleys at the sites k and I {k < i < I), with a rate equal to one, a 
layer of tiles between the two consecutive valleys is desorbed (see tile b of Fig. H]). One has 
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Figure 4: The raise and strip model. When a tile hits a peak it triggers a nonlocal 
desorption process (tile b), the tile is reflected otherwise (tile c) unless it hits a valley when 
it gets adsorbed (tile a). 



hj ^ hj -2 {j = k + l,k + 2, . . . ,1 - 1). If the tile hits a site i and 7^ 0, like in the FFM 
but unlike the RPM, the tile is reflected (see tile c of Fig. H]). 

We notice that the adsorption process is common to all three models. The differences 
are in the desorption processes. They are similar in the RSM and RPM: in general a layer 
of tiles evaporates and not a single tile. This is why one expects to flnd in the RSM model 
a value of the parameter u for which one could see conformal invariance like in the RPM. 

In the FFM and RSM the desorption takes place when a tile hits a peak and is reflected 
when the tile hits a site with no valleys or peaks. In the first model desorption is local but 
is nonlocal in the RSM. The number of active sites for desorption is equal to the number of 
active sites for adsorption. In the RPM this is not anymore the case, the number of active 
sites for desorption being the sites with Sj 7^ 0. 

The differences between the three models mentioned above will determine major differ- 
ences in their physical properties which are going to be discussed in the next sections. 

In order to study the continuous time evolution in the three models described above, 
one uses the master equation which can be interpreted as an imaginary time Schroedinger 
equation. If the system is composed by the states a = 1,2,..., Z{L), the probabilities Pa{t) 
are the solutions of the equation: 

±P,(t) = -Y.Ha,Mt). (3.1) 

The Hamiltonian H is an Z{L) x Z{L) intensity matrix: Ha^b (^ 7^ b) is non positive and 
J2aHa,b = 0. —Ha,b {cL ^ b) is the rate for the transition \b) — )■ \a). The ground-state 
wavefunction of the system |0), H\0) = 0, gives the probabilities in the stationary state: 

|0) = J2Pa\a), Pa = }}^Pait). (3.2) 

a 

In order to go from the discrete time description of the stochastic model to the continuous 
time limit, we take At = 1/(L — 1) and 

Hac = -rac (c^a), (3.3) 
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Figure 5: The average height at a site / for different values of L in the stationary state 
of the flip-flop model for u = 0.95 (under the critical point). The lattice sizes are: L = 
168,336,672,1024 and 2048. 

where r^c are the rates described above for each of the three models. 



4 The flip-flop model 

We will first discuss the stationary state of the model using some known results in combina- 
torics. We will discover in this way that we have three phases. One for u < 1, one for m = 1 
and another one for m > 1. What should we expect to find? 

If M < 1, evaporation takes over deposition of tiles. One should therefore find in the 
stationary states, in the thermodynamical limit, a finite average height. This assumption 
is confirmed by Monte Carlo simulations at n = 0.95 for various lattice sizes. The results 
are shown in Fig. [5l One can observe a constant (site independent) average height. The 
finiteness of the heights give finite values for the estimators of shared information [T3] and 
therefore we are in a gapped phase. 

For ti = 1, all the Dyck paths have the same probability and therefore can be seen as the 
paths of a restricted random walker who starts at the origin and returns after L steps. Since 
it is a random walker, this implies that one is in a gapless phase with a dynamic critical 
exponent z = 2 [T2] . 

For u > 1, Dyck paths with large heights are preferred and one can expect a growing 
interface with z = 3/2 corresponding to the KPZ universality class In this consideration 
we didn't take into account the constraints. We are going to see that our guess is not 
necessarily correct. We proceed by presenting the case m > 1 in detail. 

a) The stationary state. 

It is instructive to take L = 6 and write the Hamiltonian (13. 3p in the vector space of the 
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Figure 6: The Dyck path of the L = 6 lattice ((L + 1) points). 

five Dyck paths {Z{6) = 5). The five configurations are shown in Fig. [6l 

The configuration |1) corresponds to the substrate. Configuration |5) in which h{L/2, L) 

L/2 is the pyramid configuration. 

The matrix elements of the Hamiltonian are: 



H 



The wavefunction corresponding to the eigenvalue zero is: 
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(4.1) 



|1) + u\2) + u|3) + u^\A) + u3|5). (4.2) 

Notice that in the stationary state wavefunction each configuration gets as a coefficient 
a monomial in u with an exponent equal to the number of tiles on the top of the substrate. 
We have chosen the coefficient of the substrate (no tiles) to be equal to one. For large values 
of u, the configuration |5) corresponding to the pyramid is preferred, for small values of u, 
the configuration |1), which corresponds to the substrate has the largest probability. The 
normalization factor 

Zq{u) = 1 + 2n + + (4.3) 

is equal to the generating function for the number of tiles. One can easily check that this 
observation is valid for any number of sites and that the normalization factor Zi{u) is the 
generating function for the number of tiles for any lattice size L. Once this point is made, 
one can use results from combinatorics to get the phase diagram of the model. 

It turns out |T5] that Zl{u) is related to the Carlitz g-Catalan numbers Cn{q){L = 2n). 
The latter are defined by the recurrence relations: 

n 

Cn+iil) = E (Co = 1). (4.4) 

fc=0 

We give the first ones: 

Ci = l, C2 = l + g, C3 = l + g + 2g2 + g3^ 

C4 = 1 + g + 2g2 + 3g3 + 3g4 + 3^5 ^ (^4 5^ 
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One can define another deformation of tlie Catalan numbers 



CM = u^Cr^iu-'), (4.6) 
which are the solutions of the recurrence relations 

n 

= J2 u''Ckiu)Cr,^kiu), (Co = 1). (4.7) 

fc=0 

One has 

ZL{u) = Cniu) {L = 2n). (4.8) 
The average number of tiles is 

NL{u)=u^^^rZL{u)- (4.9) 
du 

There is no known explicit expression for Zl{u) for finite values of L and its asymptotics 
is known only for m > 1 |15] : 

]im Z2„(n) = n^/0(n-i), (4.10) 
where (j){q) is the Euler function 

oo 

0(g) = n(i-o- (4-11) 

n=l 

Using 04.91) and fl4.1Up one can compute the average number of tiles in the large L limit. 
One gets 

hm NLiu) = ^-—-^ + C{u), (4.12) 

L— >oo z 

where the L-independent term C{u) is: 

oo 7. 

Ciu) = - E (4-13) 
fc=i " ^ 

It follows that the average number of tiles is equal to those in the pyramid n{n — l)/2 for 
any u\ As a consequence, the dominant configurations are those close to the full "pyramid" 
for any u > 1. Monte Carlo simulations for finite values of L confirm this result. Taking 
u = 1.05 we show in Fig. [7^ some typical heights profiles for several lattice sizes and in 
Fig. [7)d the average values of the heights. One sees that with increasing values of L, one 
reaches the "pyramid" configuration. 

In the Appendix we consider a model in which the adsorption and desorption processes 
are the same as in the fiip-fiop model but the configuration space is not anymore Dyck paths, 
the constraint (12. 3p being relaxed. In the large L limit, the physics is the same as the one 
observed in the fiip-fiop model for u > 1. 
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Figure 7: a) Typical configurations of tlie lieigfits h{l, L) (multiplied by 2/L) in the flip-flop 
model for u = 1.05 (over the critical point) and different system sizes (L = 38, 68, 168 and 
1168). For the pyramid configuration, 2h{L/2)/L = 1. b) Average values, considering 10^ 
configurations for the same systems in (a). 

One expects that the whole u > 1 domain to be gapped and that for any initial condition, 
the systems evolve fast to reach the configurations closed to the full pyramid. We are going 
to show below that this is indeed the case. 

One can understand this result in a different way by using the Dyck paths — particle 
mapping. The desorption process peak — valley with a rate one is mapped in the hopping 
of a particle to the right: 



The adsorption process valley peak with a rate u is mapped in the hopping to the left: 



We have to take the two constraints into account. The system is confined to the size L 
and the number of particles to the left of each bond has to be larger than the number of 
vacancies. One can consider ASEP with boundaries [1] and try to mimic the constraint by 
injecting with a large rate a particles on the first site and removing particles with a rate 
P = a (the density of particles has to be equal to the number of vacancies). This picture 
should be correct ii p^ > pr {u > 1) when one has a reverse bias. The partition function 
was computed (TB] and one obtains for large systems, a factor similar to Eq. fl4.10p . 

In both cases the exponent is an area and not a length. The current also vanishes in this 
limit. There is no current in our model. 

This mapping gives absurd results if m < 1 {pr > pl)- Since we have equal densities, in 
ASEP with boundaries we are in the maximum current phase which is gapless. In our model 
there is no current, there are no density fluctuations and one is gapped. 

To sum up, the study of the stationary state suggests the following phase diagram for 
the flip-flop model. It has a gapped phase for n < 1, is gapless at n = 1 with z = 2 and 
gapped for n > 1. We now show that the time evolution results confirm this picture. 



X + 0~^0 + X, rate Pr = 1. 



(4.14) 



+ X X + 0, rate pl = u. 



(4.15) 
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Figure 8: Estimates of the dynamic critical exponent z using (14 .17^ . as a function of 1/L, 
in the three phases. 



b) Gaps and the dynamical critical exponent z 

In order to see which phases are gapped or gapless, we have diagonaUzed numerically 
(up to L = 18) the Hamiltonian (13.31) for the FFM. Since the lowest energy is zero (the 
Hamiltonian describes a stochastic process), the energy gap is given by the first excited 
state Ei{L). One expects 

Ei{L) = aiL-\ (4.16) 

where z is the dynamical critical exponent and ai a constant. A value z = implies that 
the system is gapped. In Fig. M we show the estimates for z defined as 

up to L = 18, in the three phases. One sees that zl, as expected, goes towards the value 2 
for u = 1 and to zero in the other two phases. To sum up, the phase diagram of the flip-flop 
model is as follows: 

n < 1 gapped, u = 1 gapless (z = 2), u > 1 gapped. (4-18) 

A version of the fiip-fiop model in which the rates are the same, but with a configuration 
space without the constraint (12. 3p . is presented in Appendix A. 

5 The raise and peel model - some results 

As a result of the changes occurring in the adsorption process which is nonlocal as compared 
to the FFM where it is local, the physical properties of the phases in the same phase diagram 
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also change. Like in the FFM, the u < 1 domain is gapped. Instead of a 2; = 2 phase 
transition, at n = 1 one has a z = 1 phase transition with a space-time symmetry (conformal 
invariance). For u > 1, the system stays gapless with varying dynamic critical exponents z 
{z decreases when u increases) and is not gapped like in the FFM. 

Since in the RSM we will find again a. z = 1 phase transition with conformal invariance, 
we will sum up the main results obtained in the RPM at m = 1 [7] in order to compare them 
in the next sections with those which will be seen in the RSM. 

First we discuss the stationary state. In the finite-size scaling limit l,L » 1, l/L fixed, 
the average height at a distance I from the origin, for a system of size L, has the expression: 

h{l,L) = ^\nL,, (5.1) 

where 

(5.2) 




We would like to mention that for u > 1, similar to (15. ip . h{l,L) shows a logarithmic 
increase with the size of the system L but with a different dependence on l/L. 
The density of contact points in the same limit, has the expression 

9{l,L) = ^ (5.3) 

with 

a = -^^5^HM = 0.753149.... (5.4) 

These expressions are exact. Notice the logarithmic growth of the interface and that the 
lengths dependence is all in L^. The latter is a consequence of conformal invariance. 

For M < 1, in the vicinity of the critical point u = 1, and large values of L, the average 
density of clusters K{u) is L independent and vanishes like 

if(n) = 0.596(1 - n)°-^^ (5.5) 

This result is new. It was obtained using Monte Carlo simulations. This result is a surprise. 
The scaling dimensions in the model are 1/3 and 1 [17] and no combination of them can give 
a number close to the exponent 0.78. 

The average density of peaks and valleys, in the large L limit, is 

r = l. (5.6) 

The physical meaning of (15. 6p is the following one. Adsorption takes place only on valleys 
which occupy a 3/8th fraction of the number of sites, whereas desorption takes place on a 
2/8th fraction of the number of sites. Since in a desorption process one loses more tiles 
than one gains in an adsorption process, one needs more sites with valleys than sites where 
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desorption can take place. In the RPM peaks are not active. This picture will change in the 
RSM in which sites with both valleys and peaks are active. 

The spectrum of the Hamiltonian (13. 3p . in the finite-size scaling limit, is given by 

lim = — z = 0,1,2,..., (5.7) 

L^oo ij 

where Eq = 0, Aj are the scaling dimensions and the sound velocity Vs has the value 

Vs = (5.8) 

The scaling dimensions Aj and their degeneracies {di) can be obtained from the partition 
function [17] 

oo 

i=0 n=l 

We give the first values of Aj and di. 

A = 0(l),2(l),3(l),4(2),.... (5.10) 

As we are going to see in the next Section, the RSM has a phase diagram similar to the 
RPM. For u < Uc the system is gapped, at u = Uc it is conformal invariant and gapless for 
u > Uc with varying dynamical critical exponent z. The value of Uc being not equal to one 
anymore. Which properties described above should one expect at u^? 

In the stationary state, the I and L dependence in the finite-size limit of the average 
density of contact points and heights should be again through the function Lc given by 
(15. 2p . The exponents might be different. How about the spectrum of the Hamiltonian (13.31) ? 
Since the scaling dimension Ai = 2 corresponds to the energy- momentum tensor, it should 
be not degenerate. The other values of Aj should be present but they might have other 
degeneracies. 

Using Ei{L) and (15. 7p with Ai = 2, one can determine Vg- The ratios Ei{L)/Ei{L) 
should be equal to n/L (n = 3, 4, . . .) for large L. 



6 The raise and strip model 

We are going to show that this model has a phase diagram similar to those found in the 
previous two models. A gapped phase for u < u^, a critical point at Uc and a new gapless 
phase for u > u^- In establishing the phase diagram one encounters a new problem: Uc is 
not known exactly and one has to take care of possible cross-over effects. 

We have analyzed, using Monte Carlo simulations, the heights profiles for different values 
of u and lattice sizes and found that up to u around 4.5 they are flat and the average height 
is L independent for large values of L, suggesting a gapped phase. In Fig. [9] we show for 
u = 3, the L dependence of h{L/2, L), the average height in the middle of the system. One 
sees that its value saturates fast with increasing values of L. 

The values of the average heights are larger than those observed in the RPM. For larger 
values of u there is a sharp increase of the height in the middle of the system. In Fig. [10] 
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Figure 9: The RSM. The height in the middle of the system h{L/2,L) for u = 3 as a 
function of L. The results were obtained by averaging 10^ samples. The estimated errors 
are smaller than the symbols representing the data. 

we show the heights profiles in the domain 4.5 < u < 5.4 for L = 2048. One sees that 
for u = 5.4 one sees almost a triangle with a rounded tip. A similar triangle was seen in 
the fiip-fiop model in the u > 1 domain. There is however a major difference between the 
triangles observed in the two models. Whereas in the FFM the triangle coincides with the 
pyramid configuration for which the height in the middle is equal to L/2, for u = 5.4 in the 
RSM the average height in the middle is around 30% oi L/2 = 1024. If one considers larger 
values of u the height of the triangle increases but never reaches the pyramid's height. These 
observations suggest the existence of a new phase in which the system grows linearly with 
the size of the system. We will denote this phase by LG (linear growth). We will learn more 
about it below. 

A closer inspection of the interval 4.6 < u < 4.7 suggests a phase transition between the 
gapped phase and the LG phase. In this interval the profiles are not sharp at the boundaries 
as in the gapped phase (one has an exponential fall-off in this case) and the height at the 
profile is not a triangle. One expects therefore Uc to be in this interval. 

We proceed now to a detailed analysis of the model. There is not much to say about the 
gapped phase. As we are going to show below, a good estimate for Uc is Uc = 4.685. We 
have measured using Monte Carlo simulations the density of clusters in the gapped phase 
and found for u close to Uc 



This expression is quite different of the one observed for the RPM (5.5). There are fewer 
contact points in the gapped phase in the RSM than in the RPM. This is consistent with the 
observation that the heights in the gapped phase are higher in the RSM than in the RPM. 



K{u) = 0.035(4.685 - u) 



1.733 



(6.1) 
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I/L 

Figure 10: Average height h{l, L) as a function of I in the RSM for various values of u: 4.5, 
4.6, 4,7, 4.8, 5.0 and 5.4. The size of the system is L = 2048. The average is taken over 10*^ 
independent samples. 

We first discuss the Uc physics and show that we have conformal invariance for this value 
of u. We present the results for the stationary state and time dependent observables. The 
LG phase is going to be discussed afterwards. 

a) The stationary state at u^. Conformal invariance. 

We consider the density of contact points g{l,L) for various values of u in the interval 
where we suspect to have the phase transition. If we have conformal invariance (see Section 
5), we should have in the finite-size scaling limit 

g{l,L) = CgL:^, (6.2) 

where Lc is given by (15. 2p . Cg is a constant and fi a critical exponent to be determined. The 
results of the Monte Carlo simulations are presented in Fig. [11] for L = 4096. It is shown 
that for u = 4.685 one obtains a very nice fit to the data if one takes: 

g{l,L) = 1.26L:'-^\ (6.3) 

Notice that the value of /i is very closed to 5/3. 

One can have a different look at the data for Uc = 4.685 and L = 4096 (see Fig. [T2|) 
plotting g{l, L) x L\'^^'^ as a function of sin(7r//L). One observes that within the errors, one 
obtains a constant, as expected. We have to stress that there is a small drift in the estimates 
of Uc with the size of the system. If one uses the L = 2048 data, the value Uc = 4.681 is 
preferred. This observation is illustrated in Fig. [121 

The next quantity we are looking at is the average height h{l,L). If we have conformal 
invariance, in the finite-size limit we expect to find the expression: 

h{l,L) = CHL:, (6.4) 
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ln(Lsin(Til/L)K) 

Figure 11: In [g{l, L)] in the RSM as a function of Lc for L = 4096 and for the values of: 
u = 4.681, 4.682, 4.689. The results are obtained by averaging over 10^ Monte Carlo steps. 
The fitting curve ai u = 4.685 is y = —0.234 — 1.65x. 
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sin(Kl/L) 

Figure 12: The density of contact points g{x,L) multiplied by (Lc)^'^^ as a function of 
sin{7Tl/L) for u = 4.685 and L = 4096. We show the results of five samples with 10^ Monte 
Carlo steps 
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Figure 13: Effect of the lattice size on the determination of Uc for the lattice size L = 2048 
by averaging 10^ Monte Carlo steps. The value Uc = 4.681 gives a better estimate for Uc 
than Uc = 4.685 obtained for L = 4096. 

where z/ is a critical exponent and Ch is a constant. Using the data obtained for Uc = 4.685 
in Fig. [14] we plot L) x Lj''^ as a function of sin(7r//L) and find almost a constant value. 
This implies u = 0.5. We have measured u considering h{L/2, L) for various values of L and 
our best estimates are 0.50 < u < 0.52. An exact estimate for u is hard to get since, as 
mentioned above, the estimates of Uc change slightly with L. 

Finally one has measured the density of peaks and valleys for large values of L and 
obtained 

T ^ 0.35. (6.5) 

What have we learned up to now about the Uc phase transition? The data are compatible 
with conformal invariance. We will also show below that the critical exponent z = 1 and 
that the spectrum of the Hamiltonian is what we expect. Although they are both conformal 
invariant, the phase transitions at u = 1 for the RPM and at Uc for the RSM are different. 
The average height growth logarithmically in the RPM and like a power in the RSM. The 
density of contact points has an exponent 1/3 in the RPM and probably 5/3 in the RSM. 
Also the density of peaks and valleys are different. There are fewer peaks and valleys in 
the RSM. This can be understood as follows: since the phase transition takes place at a 
large value of u, in the stationary state, many tiles are adsorbed therefore many have to be 
desorbed. This implies that one needs many sites without peaks or valleys. 

b) Time dependent phenomena at Uc- 

In order to determine the value of the dynamic critical exponent we use the properties 
of the Family- Vicsek scaling function [18]. We first consider the density of clusters K{L,t). 
From (6.3) follows that the number of clusters is finite therefore the density of clusters in 
the stationary state behaves like 

K{L) = Ck/L, (6.6) 
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Figure 14: The average height h{l,L) multiphed by LJ^^'^ as a function of sin(7r//L) for 
u = 4.685 and L = 4096. The resuhs shown are obtained by averaging 10^ Monte Carlo 
steps. The estimated error are smaller than the symbols in the figure. 

where Ck is a constant. This implies that for large values of L and short times one expects 

m = % (6.7) 

being Dj^ a constant. In Fig. [15] we show the short time dependence of K{t\ and a fit to 
the data gives 

\^\K{t)\ = -1.38 - 0.99 ln(t), (6.8) 

from which we get the value z = 1.01 extremely close to the value z = 1 required by 
conformal invariance. 

A similar analysis was done for the average height at the half value of L, h{L/2, L). Since 
in the stationary state h{L/2,L) ^ L^/^, we expect, if z = 1 that for short times and large 
values of L to find h{L/2, L,t) ^ t^/^. This expectation is confirmed by the data shown in 
Fig. [ini where for short times the fit gives: 

ln[/i(L/2, L, t)] = 0.154 + 0.49 ln(t). (6.9) 



A final confirmation of conformal invariance at Uc is obtained looking at the first excited 
states of the Hamiltonian (13. 3p at Uc = 4.685. We have diagonalized numerically the Hamil- 
tonian up to L = 18. In Fig. [T7]we show for convenience L x Ei{L)/2'k and not L x Ei{L)/'k 
as a function of 1/L. Using the data and (5.7) with A = 2 one gets Vs ~ 3.7 (a value of 
the same order as the value (15.81) in the RPM). We also show the ratios E2{L)/Ei{L) and 
E2,{L) / Ei{L) , as a function of L. The data nicely converge to the value 3/2 indicating that 
we have two levels with A = 3 (in the RPM one has only one level with this value, see 
(5.10)). 
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Figure 15: Determination of the dynamic critical exponent z at u = 4.685. The density 
of clusters versus time on a log scale. The substrate was used as the initial condition, and 
the average was taken over 10^ samples. The fitting y = —1.38 — 0.99a; was obtained in the 
region inside the ^ ^ 
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Figure 16: The average height at L/2 as a function of time for L = 512, 1024, 2048, 4096 on 
log scales. The initial state was the substrate, and the average was over 10^ samples. The 
region where the fitting was taken is also shown. 
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Figure 17: Finite-size scaling analysis of the spectrum of the Hamiltonian at Uc = 4.685. The 
data are obtained by diagonalizing the Hamiltonian up to 18 sites. LEi{L)/2'it, E2{L) / Ei{L) 
and E'i{L) / Ei{L) are shown as a function of 1/L. The extrapolated values for Ei/{2tiL) is 
also shown. 

The existence of two levels with A = 3 is a puzzle. Since the central charge of the Virasoro 
algebra is c = (for a stochastic model the ground state of the Hamiltonian is always equal 
to zero), the character corresponding to the vacuum representation is Xo{q) = ^{q) (see 
(15. 9p ). There is only one scaling dimension with A = 3 (see (I5.10p ). On the other hand, 
A = 3 is not in the Kac table according to which the scaling dimensions are 

A^,, = [(3r-2s)2-l]/24 (6.10) 

with r and s positive integers and therefore one has to add to the spectrum a standard 
Virasoro character: 

oo 

x3(g)=g'n(i-o-'- (6.11) 

71=1 

In the finite-size scaling limit, the spectrum of the Hamiltonian is given by the partition 
function Xo{^) + Xs{q)- This implies the following values of Aj and degeneracies {di): 

A = 0(1), 2(1), 3(2), 4(3), 5(4), .... (6.12) 

In order to check if the levels with the corresponding degeneracies are seen in our model, 
we have looked at the lower part of the spectrum of the Hamiltonian for L = 18. Normalizing 
the eigenvalues by taking £'i(18) = 2, one gets the following values: 

0; 2; 3.16; 3.24; 4.35; 4.36; 4.54; 5.20; 5.4; 5.6; 6.0; .. . (6.13) 

We have checked that the values of Ei{L) / Ei{L) converge from above from larger values and 
cluster towards their asymptotic ones. We conclude that, although the lattice sizes are up 
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Figure 18: The u dependence of H{u). Monte Carlo simulations for lattice L = 16384. The 
results were obtained by taking 10^ Monte Carlo steps. The estimated errors are smaller 
than the symbols. 

to L = 18, at least for the lower part of the spectrum there is agreement between f l6.13p and 

We have tried to see if at Uc, the Hamiltonian has hidden symmetries for finite L. They 
could show up in the observation of degeneracies of the spectra at finite L. In fl6.13p one 
notices that two couples of values are very close (3.16 and 3.24 respectively 4.35 and 4.36). 
Since Uc is not known exactly we have checked if small changes in the value of u around the 
value 4.685 could make the levels degenerate for all values of L. Our investigation gave a 
negative result. Degeneracies were observed only for one of the two couples and only for one 
value of L. This observation probably excludes the existence of Jordan cells in the finite-size 
scaling limit since in all known examples the Jordan cells occurring in the conformal theory 
appear already for finite values of L. 

c) The linear growth phase (u > U(.). 

As mentioned earlier, in this phase the heights profiles in the stationary state are triangles 
(see Fig. [10] for L = 2048) with heights which are fractions of the maximum possible height 
L/2. The maximum height is seen in the stationary states of the FFM. We have studied the 
u dependence of 

H{u) = 2h{L/2,L)/L (6.14) 

using Monte Carlo simulations for a very large lattice (L = 16384). The results are shown 
in Fig. [m One notices, as expected, that H{u) increases with u (for u infinite, H{u) = 1, 
the triangle becomes the pyramid). 

Since H{u) has a clear geometrical meaning, it can be seen as an order parameter in the 
LG phase [I9]. 

In order to determine the values of the dynamic critical exponent z, we have performed an 
analysis similar to the one done for u = Uc (see eqs. (16. 6 p and (16. 7p ). We have first convinced 
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Figure 19: The number of clusters LK{L,t) as a function of time for two different initial 
conditions (pyramid in red and substrate in black). The lattice size is L = 4096 and u = 5. 
On each evolution there were 4 x lO^samples in the Monte Carlo simulations. 



ourselves that for m > 1 in the stationary state and the large L limit the number of clusters 
stays finite. Starting with two different configurations (the substrate and the pyramid), we 
have looked at the time dependence of the number of clusters Cxif) and checked if they 
converge at the same value Ck- The results of the computer simulations for L = 4096 and 
u = 5 shown in Fig. [T9] show that this is indeed the case. This observation implies that in the 
stationary state, the density of clusters is given by eq. (16. 6p . The next step is to determine 
z using eq. (16. 7p . 

In Fig. |20]we show the short time dependence of the density of clusters for various lattice 
sizes and u = 7. One has used the substrate as initial condition. Notice the data collapse 
which allows to determine z. A similar analysis was done for other values of u. Our estimates 
for z are 0.77 {u = 5.5), 0.66 {u = 6), 0.5 {u = 7), 0.43 {u = 8), 0.35 {u = 10). 

We conclude that in the LG phase one is gapless, the dynamic critical exponents decrease 
in value with increasing values of u. We had a similar behavior in the RPM. 



7 Conclusions 

We have considered three one-parameter dependent stochastic models (FFM, RPM and 
RSM) defined on Dyck paths. Dyck paths have fixed boundaries (12. 2p and are confined to 
the upper half plane (12. 3p . The parameter denoted by u in the text, is equal to the ratio of 
the adsorption and desorption rates. Adsorption is local in the three models, desorption is 
local in the FFM and nonlocal in the other two. The phase diagrams of the three models 
are similar: for u < u^. the model is gapped, it undergoes a phase transition at u^. where it 
is gapless, for n > it is again gapped for the FFM but gapless for the RPM and RSM in 
which desorption is nonlocal. 
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Figure 20: Short time dependence of the density of clusters for several lattice sizes of the 
RSM at u = 7. It is also shown the fitting results in the region of the dashed rectangle for 
L = 2048 and L = 4096. On each evolution there were 4 x 10^ samples in the Monte Carlo 
simulation. 



The aim of the paper was to understand the effect of having nonlocal rates. The raise and 
peel model was shortly described here because it is, as far as we know, the first intensively 
studied model [5,6] in which the desorption process, which takes place when a tile hits a 
slope, is nonlocal. The peaks are not active. In the raise and strip model presented for the 
first time in this paper, the slopes are not active but when a tile hits a peak it triggers a 
nonlocal desorption process. In order to clarify the role of nonlocality in the latter model, 
we have presented the flip-flop model which is similar to the RSM (the peaks and valleys are 
the active sites) with the major difference that desorption is local. 

Two Dyck paths play an important role in these models: the substrate (no tiles) and 
the pyramid configuration in which the height at the middle is equal to half the size of the 
system. If n = 0, the stationary state is the substrate configuration in all three models. 

In the flip-flop model Uc = 1 and the dynamical critical exponent z = 2 (one has a random 
walker). We have studied the spectra of the Hamiltonian for u < 1 and n > 1 and obtained 
that the system is gapped in both cases. In the stationary state, the average height is finite 
for u < 1, and increases like L^/^ at u = 1 (L is the system size). For any m > 1, in the large 
L limit, the system has small fluctuations around the pyramid configuration (the average 
height at the middle of the system h{L/2, L) is equal to L/2). This result was obtained using 
known results from combinatorics. A variant of this model in which the configuration space 
is changed (RSOS paths without the restriction (2.3)) and the rates are local is presented in 
Appendix A. The results are similar to those obtained for the flip-flop model. 

In the raise and peel model, Uc = 1, the dynamical critical exponent z = 1 and one has 
conformal invariance. It is this special property which inspired our work presented here. Is 
conformal invariance a consequence of nonlocality? That this might be the case is suggested 
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by the observation that in the peak adjusted raise and peel model [TT] where adsorption is 
also nonlocal, one has conformal invariance. In the stationary states the average heights are 
finite for u < 1 and increase logarithmically with L not only for u = 1 but also for m > 1. 
There are many properties of this model which are known exactly because the Hamiltonian 
is integrable at u = Uc = I- 

The raise and strip model is most probably not integrable. All our results are based 
on Monte Carlo simulations on large lattices. The best estimate for the critical point is 
Uc = 4.685. For u < Uc and large L, the average values of the heights stay finite albeit 
larger than in the RPM. For u = Uc, one obtains z = 0.99 very close to the value z = 1. 
Several tests suggest that one has conformal invariance. At Uc the RSM model has different 
properties than the RPM. The finite-size scaling spectrum of the Hamiltonian are given by 
different Virasoro modules in the two models. The critical exponents are also different. The 
average height increases like a power of L and not like In L. For u > Uc the dynamical critical 
exponent z decreases if u increases, similar to what was observed in the RPM. The analogy 
stops here however. Whereas in the RPM the average height increases logarithmically with 
L, in the RSM the heights profile is very different. For large values of L, the profile is 
a triangle with the tip increasing linearly with L. We have therefore called this domain 
of M, the "linear growth phase". For u very large the triangle gets close to the pyramid 
configuration seen in the fiip-fiop model. 

We believe that the phase diagram observed in the RPM and RSM is general for the 
class of models of nonlocal growth. This has of course to be proven looking at other models. 
Such a study might also bring a better understanding of some properties of the RSM (the 
values of the exponents are just one example). 

The relevance of our work for physical systems like polymers with nonlocal interactions 
[m |22l [23] remains to be studied. 
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A Filling a square with tiles 

This is the fiip-fiop model described in Section 3 in the configuration space of RSOS paths 
without the restriction f l2.3p . The paths are allowed to move in the lower half-plane and it 
is amusing to see what is the effect of changing the configuration space. If m = 1, one has 
a random walker with fixed ends at positions and L. It is convenient to see the paths as 
describing an interface between tiles which fill part of a square and a rarefied gas of tiles. In 
Figure [21] we show such a path in the case L = 6. Out of a maximum of 9 tiles, the square 
is filled with 7 tiles. 

To see how the model works, we take L = 4. There are 6 configurations in this case as 
shown in Figure [25] 
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Figure 21: Configuration space for L = 6. Tlie interface is a path which separates the seven 
tiles inside the square from a rarefied gas of tiles. 




Il> I2> I3> |4> I5> I6> 

Figure 22: The six configurations for L = 4. 



The Hamiltonian obtained using the rules of the fiip-fiop model is 
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(A.l) 



The wave function corresponding to the eigenvalue zero is: 



n^|l) + n3|2)+n2(|3) + |4)) + n|5) + |6). 



(A.2) 



In the stationary state, each configuration gets as a coefficient a monomial in u with an 
exponent equal to the number of tiles inside the square. The partition function is 



Za{u) 



1 + u + 2n^ + + u^. 



(A.3) 



One could have obtained directly flA.2p using the matrix-product method [20] for an open 
system with the processes given by (4.14) and (4.15). If a step up in the path is given by a 
generator U and a step down by a generator D, we can use the algebra 



UD = uDU 



(A.4) 



to find for example, that the coefficient of the configuration |1) [UUDD) gets a factor 
compared with the configuration |6) [DDUU). 
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(A.5) 
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(A.6) 



It is easy to show that the number of tiles inside the square 



Nl{u) = u—\iiZl{u) 



(A.7) 



in the large L limit is: 



ri^ + C{u) if M > 1 
nV2 if u = 1 

-C{l/u) ifn<l 



(A.8) 



where C{u) is given by f l4.13p . This result is amusing since it shows that for any m > 1, in 
the large L limit, the important configurations are the same as for the fiip-fiop model. For 
M < 1, the relevant configurations are those near the empty square. They are obtained by a 
simple reflexion (top — )■ bottom) of the configurations relevant for m > 1 (the Hamiltonian 
is invariant under the transformation n — )■ 1/u with a change of the time scale). At u = 1, 
one has the random walker. 

The phase diagram of this model is the same as the one of the fiip-fiop model. The only 
difference between the two models in the large L limit is the heights profiles in the u < 1 
domain. 
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